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and the random first-order transition theory 

Laura Foini/ Florent Krzakala,^ and Francesco Zamponi^ 

^ Universite Pierre et Marie Curie - Paris VI, LPTHE, 

Tour 13 Seme etage, 4 Place Jussieu, 75252 Paris Cedex 05, France 

^ESPCI ParisTech, CNRS UMR 7083 Gulliver, 10 rue Vauquelm, 75005 Pans, France 

'^LPT, Ecole Normale Superieure, CNRS UMR 8549, 24 Rue Lhomond, 75005 France 

In this paper we revisit and extend the mapping between two apparently different classes of 
models. The first class contains the prototypical models described -at the mean-field level- by 
the Random First Order Transition (RFOT) theory of the glass transition, called either "random 
XORSAT problem" (in the information theory community) or "diluted p-spin model" (in the spin 
glass community), undergoing a single-spin flip Glauber dynamics. The models in the second class 
are Kinetically Constrained Models (KCM): their Hamiltonian is that of independent spins in a 
constant magnetic fleld, hence their thermodynamics is completely trivial, but the dynamics is such 
that only groups of spin can flip together, thus implementing a kinetic constraint that induces a 
non-trivial dynamical behavior. A mapping between some representatives of these two classes has 
been known for long. Here we formally prove this mapping at the level of the master equation, and 
we apply it to the particular case of Bethe lattice models. This allows us to show that a RFOT 
model can be mapped exactly into a KCM. However, the natural order parameter for the RFOT 
model, namely the spin overlap, is mapped into a very complicated non-local function in the KCM. 
Therefore, if one were to study the KCM without knowing of the mapping onto the RFOT model, 
one would guess that its physics is quite different from the RFOT one. Our results instead suggest 
that these two apparently different descriptions of the glass transition are, at least in some case, 
closely related. 



I. INTRODUCTION 

The Random First-Order Transition (RFOT) ^^ and 
the dynamic facilitation theory based on Kinetically Con- 
strained Models (KCM) [3-01 are often viewed as alter- 
native pictures of glass formation [6l. [8l-[lO|. Indeed, in 
the former the glass transition is induced by an underly- 
ing thermodynamic transition, the so-called "Kauzmann 
transition" , while in the latter the thermodynamics is as- 
sumed to be completely trivial, and the glass transition 
is a purely dynamical arrest induced by complex kinetic 
constraints^. 

The RFOT theory^ is based on the analysis of a par- 
ticular class of mean-field spin glass models with p-spin 
interactions [l[ . It has been known for long [5|, [Il|, [ij] 
that some particular finite-dimensional version of these 
p-spin models (the so-called "plaquette models" ) can be 
explicitly mapped onto KCMs. However, the physics of 
these models, at least in two dimensions, seems to be 



Note that sometimes (especially in the mathematical literature) 
a more stringent definition of KCM is employed. Namely, one 
requires that the dynamics of a finite size system has absorbing 
(or "blocked") states. In the following, we will just define a KCM 
as a system that has a trivial thermodynamics, and a complex 
dynamic behavior induced by specific dynamic rules. 

In this paper, by RFOT theory we denote the ensemble of the- 
ories that can be related to the exact solution of mean-field spin 
glass models (therefore including mode-coupling theory, molec- 
ular liquid replica theory, and finite-dimensional extensions). In 
particular, we will mostly be concerned with mode-coupling-like 
dynamics, but we will keep calling it RFOT dynamics. 



quite different from the RFOT one d [TT|-[l3|. thereby 
preventing one to establish a direct mapping between the 
two approaches. 

The aim of this work is to discuss an example in which, 
on the contrary, a KCM can be exactly mapped into a 
model described by the RFOT theory. In order to do 
this, we first generalize the mapping discussed in [5|, llll - 
|l3| to an arbitrary geometry, with the only requirement 
that i) each spin is involved in exactly p interactions, 
a) each interaction involves exactly p spins, and Hi) the 
ground state is unique. By using the method developed 
by Cosset [lj| for quantum Hamiltonians, we show that 
the Fokker-Plank operator of the p-spin model can be 
mapped into the one of the KCM. 

We then focus on a Bethe lattice geometry, because 
this case allows for a very complete and detailed anal- 
ysis, and in particular the Bethe lattice p-spin model is 
exactly described by the RFOT theory. More precisely, 
we consider the two following models: 

• The first model is a spin model defined on a random 
regular factor graph (or Bethe lattice) geometry, 
where each spin is randomly connected to three in- 
teractions and each interaction involves three spins. 
The dynamics is a single spin-flip Clauber dynam- 
ics. The model is known as "random 3- XORSAT" 
in the information theory community, and as the 
diluted 3-spin in the spin glass community. It 
was studied intensively by both communities |15l - 
Il9| , and it was shown that the dynamical spin- 
correlations show a plateau in time that is natu- 
rally interpreted from a thermodynamical point of 
view in terms of the emergence of disconnected er- 



godic states that trap the system's dynamics (see 
Fig. [1}. Rigorous resuhs have in fact shown that 
for such a system the Gibbs-Boltzmann measure 
can be decomposed into disconnected components 
and purely dynamical quantities are recovered from 
a thermodynamical study [13, H, [20| • We will refer 
to this model simply as "XORSAT" in the follow- 
ing. 

• The second model is defined on the same geome- 
try as the first one. Its Hamiltonian is given by 
non-interacting spins in a magnetic field, and its 
dynamics only allows a simultaneous flipping of a 
group of three spins connected to the same "interac- 
tion" . This model displays the typical behavior of 
a Bethe lattice KCM 21], i.e. it has trivial equilib- 
rium thermodynamical properties but the dynamic 
rules lead to a complex dynamical behavior char- 
acterized by a finite-temperature dynamical tran- 
sition. As in other KCMs, despite the simplicity 
of the equilibrium Gibbs-Boltzmann measure, the 
dynamical correlations display slow relaxation and 
characteristics that resemble those of glassy sys- 
tems. In the following we will refer to this second 
model as the KCM. 

Thanks to the general mapping discussed above, we show 
that these two models are in fact completely equivalent: 
the master equations describing their dynamics can be 
mapped one into the other by a suitable transformation, 
that maps the groups of three interacting spins ( "plaque- 
ttes" ) of XORSAT into the spins of the KCM. We show 
that while different observables share a similar dynamical 
behavior, and in particular a plateau in their time cor- 
relations associated with the dynamical transition, they 
crucially differ once that they are chosen as order param- 
eters in a static computation. 

Therefore, some observables fail to identify the ergod- 
icity breaking. Remarkably, this is the case for all local 
observables in the KCM: detecting ergodicity breaking in 
the KCM therefore requires the use of non-local (hence 
very unnatural and complicated) observables. To illus- 
trate this point, we will consider two distinct order pa- 
rameters for the XORSAT model: the spin overlap Qs 
and the plaquette overlap qp. We will show that qs is 
a good order parameter, while qp is not. The plaque- 
tte overlap for XORSAT becomes the spin overlap for 
the KCM, therefore the spin overlap is not a good order 
parameter for the KCM. Moreover, the spin overlap for 
XORSAT becomes a very complicated non-local function 
for the KCM; therefore, detecting the glass transition by 
means of a static computation in the KCM requires the 
use of a very unintuitive order parameter. If one were to 
study the KCM without knowing of the underlying map- 
ping on XORSAT, it would be impossible to identify the 
correct order parameter, and one would conclude that 
the dynamic transition cannot be detected via a static 
computation, and that the dynamical facilitation picture 
is completely different from the RFOT one. 



Our results suggest that this conclusion might be 
wrong: at least for the model we investigated, the two 
scenarios are dual to each other. The same could be true 
in more general models 22] . 

The paper is organized as follows. In section |TT] we 
establish in general the mapping between the Fokker- 
Planck equations describing the two classes of models. 
In section Hill we focus on the Bethe lattice models and 
present the results for both dynamic correlation functions 
and the thermodynamic RFOT potential. Finally, in sec- 
tion |IV] we discuss our results and their implications. 



II. MAPPING p-SPIN MODELS INTO KCMS 

In this section we show how to characterize the equi- 
librium dynamics of the models and how to construct the 
mapping between the two. 

Let us first outline the general Fokker-Planck equation 
for the stochastic dynamics of N Ising spins, Ui = ±1 
and z = 1, . . . , iV, described by a generic energy function 
E{a_). We define P{q_,t) the probability that at time t 
the system is in the configuration a. Suppose that the 
elementary dynamics is induced by the flipping of a set 



of variables k = {4,*fe,''' , *fe } of size |fc|, with rate 
^is. iS): where a}^' represents the configuration that 
is obtained from a by the flip of the set k. Then, the 
probability P{g_, t) satisfies the equation: 



dt 



Y.w{a,s.^-^))P{q}^\t), 



(1) 



where the sum over k stems for all possible sets of vari- 
ables. In order to reach equilibrium the rates have to 
satisfy the detailed balance condition: 



w{^^^\^)e 



-lSE{a) 



w{a,a 



W^p-/3S(2'"') 



(2) 



and in the following we will adopt the common choice 



z/;(a«,a)=u;oe-^(^(^*'')^^(^» 



(3) 



with Wo = 1 in order to set the unit of time. 

Defining the matrix H such that {o_\H\a') — 
Efe«^(a«,a)5(a,a') - Efc «^fc^^'^)'^(^^'^^'), and 
\P{t)) the probability vector such that (a;|P(i)) — 
P{g.,t), Eq. dl]) can be written as a Schrodinger equa- 
tion in imaginary time: 



dlPjt)) 
dt 



-H\P{t)) 



(4) 



It is often convenient to introduce the transformation 
induced by the diagonal matrix [/, such that {q\U\a!) — 



e2^(2:)(5(CT,CT'), which turns Eq. (g]) into: 



d\P{t)) 
dt 



-H\P{t)) 



(5) 



with \P{t)) = U\P{t)) and H = UHU'^ The advantage 
of this formulation is that the matrix H is Hermitian. 

Under the condition that the rates satisfy Eq. Q, one 
can verify that the equihbrium Gibbs-Boltzmann mea- 
sure P{a,t) — Pgb{^) oc exp{—l3E{a)) is a stationary 
solution of Eq. ([T]) and it represents the ground state of 
the Fokker-Planck operator H in Eq. ((4]), and the same 
holds for Pgb{2.) oc exp{-PE{a)/2) and H. 



B. Kinetically constrained model 

We now introduce the kinetically constrained model 
(KCM), using different notations that will be suitable 
for the mapping to be discussed in the next section. We 
consider a system of N Ising spins Sq = ±1 whose energy 
function reads: 



N 



E 



KCM 



u) = -E^ 



(8) 



A. XORSAT model with spin-flip dynamics 

We consider a model of N Ising spins characterized by 
3-body interactions with Hamiltonian: 
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(6) 



where a labels the interactions, each one involving a 
triplet of spins. Such Hamiltonian represents a partic- 
ular case {p = 3) of a general class of spin models with p- 
body interactions that in computer science is also known 
as XORSAT problem. 

The graphical representation of the problem can be ob- 
tained through the introduction of a "factor graph" , by 
assigning a "variable node" to each spin i and an "inter- 
action node" to each interaction a and connecting every 
spin to the interaction nodes in which it takes part. Here 
we assume that each spin is connected to exactly three 
interactions, hence the connectivity of spin variables is 
c = 3. For the Hamiltonian in Eq. ©, the connectiv- 
ity of each interaction node is p = 3 and the number M 
of these nodes must satisfy Mp — Nc, then M — N. 
In particular, from Eq. (j6|), the ferromagnetic configura- 
tion ai = 1 Wi = 1, . . . , N minimizes all the interactions 
a — 1, ... ,M and in what follows we will restrict to the 
case where this is the only ground state (we will come 
back to this point in section lin|) . 

For this model we study the dynamics induced by sin- 
gle spin flips, so the sum over k in Eq. ([T]) is simply 
the sum over all the spins and the rates w{a}^',a) con- 
nect configurations differing only for a single spin ctj. 
Then, according to Eqs. ^ and ([6]), wxoRi2.^^\^) = 
e '^^"^^''^•i'^'a'^'a J where here and in the following a G 
di indicates the neighborhood of i, i.e. the interactions 
a to which the site i is connected. Performing the opera- 
tions outlined in the previous section, one arrives to the 
following form for the Hermitian Fokker-Planck operator 
that appears in Eq. ([5]): 



H 



N 

XOR — 2^ ^ 
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(7) 



Due to the non- interacting nature of Ekcm {§) the ther- 
modynamical properties of the system are independent 
of the geometry where the system is defined. 

In order to specify the constrained dynamics, we as- 
sume that the model is defined on a regular random 
factor graph, with p = c — i, exactly as before, and 
we define the dynamics by the flip of triplets of spins 
{a\,al,a'l} surrounding the same "interaction" node i 
(hence we use a dual notation with respect to the XOR- 
SAT model). The set k in the Fokker-Planck Eq. ^ is 
now the set of the three spin identified by i and, follow- 
ing our choice in Eq. Q, the rates read wkcm{s^^\s) = 
exp(— /3^^g^j Sa). The Hermitian evolution operator 
that enters in the Fokker-Planck Eq. ([5]) for the KCM 
can be expressed in terms of Pauli matrices as follows: 



H 



KCM 
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(9) 



Note that the choice of the indices z and x in Eq. ^ is 
just an arbitrary choice of basis for the representation of 
the spin configurations and one can safely exchange the 
two indices. 



C. Mapping 

Following the procedure outlined in |Tj] in the context 
of quantum spin systems, in this section we describe the 
mapping between Hxor, defined on a regular graph with 
c ^ p and such that Exor has a unique ground state, 
into Hkcm defined on the dual graph of the interaction 
nodes. 

Consider the operator Hxor. defined in Eq. ([7]) for a 
given graph with a unique (ferromagnetic) ground state. 
Thanks to this property, for each interaction a, it exists 
one and only one configuration of spins ct° minimizing all 
the interactions but a. To show this, let us introduce the 
N X N adjacency matrix A such that Aai = 1 if the spin 
ai enters in the interaction a and Aai — otherwise, and 

-. The matrix A is invert- 



where df'^'^ are the standard Pauli matrices. 



define y°' such that yf = 
ible because by assumption the system Ay = (that is 
equivalent to the condition of minimizing Exor{S-)) ^^^ 
a unique solution, corresponding to y = 0. Defining the 
vector e° such that e^ = 



one gets y"^ — A e" 



Then, one can introduce the following operators, de- 
fined on the interaction nodes of the underlying graph of 
Eq. 0: 



-^a ^ i'- i'^ i"^ ? 



(10) 



These are well defined spins operators because they sat- 
isfy the correct commutation relations, [s^jS^] = if 
a y^ b, {slJl} = and (s^)^ = (s^)^ = I, where I is 
the 2x2 identity matrix. Moreover, it turns out that 
defining e* such that e* = Sij, and using the relation 






'Y.st = Y.f (11) 



the mapping above can be inverted as follows: 

N 



(12) 



Then, enforcing the mapping in Eq. (jlOp and (|T2j) in 
HxoB, one obtains the operator Hkcm of Eq. ([9]) de- 
fined on the dual lattice. Note that here we discussed the 
mapping for p = c = 3, but the extension to any other 
value of p = c is straightforward [14]. 

We have therefore shown that HxoR defined on a given 
graph can be written as Hkcm defined on the dual graph 
by means of the transformation p2|) . and viceversa. Note 
that this implies in particular that the ground state of 

HxOR, 
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(13) 



whose square is the Gibbs distribution, can be written as 
the one of Hkcm on the dual graph: 



{s.\GSkcm) 



1 



g-|/3£KCAf(£) 



(14) 



and moreover the partition functions of the two mod- 
els (that correpond to the normalization of the ground 
states) are identical: 



Z 
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-0Ekcm{s) 



(15) 



graphs characterized by N spins and M interactions such 
that each spin appears exactly in c = 3 interactions and 
each interaction connects precisely p = 3 spins, and such 
that Eq. ^ has a unique ground state (the last condi- 
tion is needed in order to use the mapping discussed in 
section IIICp . We will be interested in the properties of 
the model in the thermodynamic limit N = M ^ oo. In 
this limit, such graphs have the property of being locally 
tree-like, as a Bethe lattice should: indeed, typical loops 
have length of order log N. The lack of a finite dimen- 
sional geometry make statistical models defined on such 
lattices intrinsincally mean-field like. Most importantly 
this property allows for an analytical treatment via the 
so-called "cavity method", and many results have been 
derived in this way |17| . The thermodynamical prop- 
erties of Eq. ^ are self-averaging with respect to the 
ensemble of the graphs so defined and it turns out that 
in the case considered here, k — c, the ground state is 
non-degenerate with finite probability j23j. Hence, the 
restriction to graphs that have a unique ground state, 
that has to be made here in order to derive the map- 
ping, does not alter the main properties of the model. 
Note that the ensemble of regular random graphs with 
c = p is self-dual, hence the KCM will be defined on 
the same ensemble of graphs. Note also that the two di- 
mensional triangular lattice considered in '5'| represents 
a (very atypical) instance of such ensemble of graphs. 



A. XORSAT model 

The XORSAT model on the random regular factor 
graph with p — c — 3 follows exactly the RFOT scenario 
of the glass transition [13, H, [IJ • It is characterized by 
a dynamic transition at a temperature Td — 0.5098(2), 
below which ergodicity is broken and the equilibrium dy- 
namical correlation functions of local observables never 
relax. In the RFOT scenario there is a second tempera- 
ture Tfe at which a thermodynamic transition to a replica- 
symmetry broken state takes place 'l'|. However, here 
Tk = 0, because of the mapping onto independent spins 
that ensures that the free energy is analytic at all tem- 
peratures, see Eq. ([T5|): this can also be checked by an 
explicit replica or cavity computation [15nl7. 24 , 25] . 



1. Spin correlations 

In the formalism used above the local spin-spin corre- 
lation function for the XORSAT problem reads: 



III. BETHE LATTICE MODELS 

Here we specialize our discussion to Bethe lattice mod- 
els. Different definitions of Bethe lattice have been used 
in the literature, and here we consider as a model of Bethe 
lattice a random regular factor graph: we construct an en- 
semble of graphs by giving uniform probability to all the 



1 ^ 
C^f''''(i) = i^E('^*W^»(0))' 



^tt (16) 

where {z.\GSxor) is the ground state of Hxor defined 
in Eq. (|13p . These correlation functions are shown in 
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FIG. 1: Local spin-spin correlation function C^ (t) for the 
XORSAT model at different temperatures, from left to right 
r=0.7, 0.65, 0.6, 0.58, 0.56, 0.55, 0.54, 0.535, 0.53, 0.525, 
0.52, 0.515, 0.51, 0.505, 0.5, 0.48, 0.45. The inset shows the 
comparison of the long time dynamics with the spin overlap 
that is obtained through the cavity method, which is reported 
as a dashed line for each temperature [13, [23] . 
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Fig. [TJ ioT T < Td, they tend to a plateau q\ 
long times. The dynamical correlation functions, here 
and in the following, are obtained using the Metropolis 
algorithm for a system of iV = 60000 spins. Sampling 
the initial configuration below Td is extremely hard with 
standard techniques, but it can be easily achieved using 
the planting technique [13, HI 113 • 

The plateau value qf~'^ can be recovered from a sta- 
tistical treatment of the thermodynamical properties of 
Eq. ([6]) via the cavity method 17, 25]. In this treatment, 
it is assumed that the system at long times decorrelates 
inside one of the typical ergodic components a at tem- 
perature T; then 
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(17) 



where the overline indicates the equilibrium average over 
ergodic components. The equivalence between these two 
quantities is shown in the inset of Fig. [1] 

A more detailed correspondence between statics and 
dynamics can be established by means of the Franz-Parisi 
thermodynamical potential [23,1231, defined as: 
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(18) 



overlap q from the spin system defined by a_. For mean- 
field models described by the RFOT theory, the potential 
V^^^{q) has the following properties, see Fig. [2] 

• It is a convex function with a single minimum at 
g = for T > Trf 



• For Td > T > Tk it displays a secondary minimum 
at gf °^ > 0, such that ^/^^(gf ^fl) > y/ofl,(o). 
The difference between the two minima can be in- 
terpreted as the free energy cost to force the two 
system to have a large overlap. 

• For T < Tk the secondary minimum becomes ther- 
modynamically stable, V.^^^iq^^'^) < y/^-^(0), 
signaling the onset of replica symmetry breaking. 

For the model under consideration, i.e. XORSAT with 
p = c, it turns out that T^ = 0, as we already discussed; 
thus the third regime is not relevant for the present dis- 
cussion. Importantly enough, the value of the overlap 
q^OR ^i^QYe the secondary minimum is attained when 
T < Td represents the plateau at which the dynamical 
correlations of the spins converge in the long time limit. 
Thus, the non-trivial dynamical behavior can be under- 
stood in terms of thermodynamical quantities, which is 
the key feature of RFOT theory [H-ls^. 



2. Plaquette overlap 

Let us now turn to the discussion of a second observ- 
able for the same model. Instead of focusing on spin-spin 
correlations we study here the dynamics of a "plaquette" 
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where ^^(ct, ct') — -^ X]i=i '^i'^'i i^ t-^i^ spin overlap func- 
tion that measures the similarity of two configurations, 
(J and <j' . The potential V^'~'^{(i) represents the free en- 
ergy of a system of spins a;' constrained to be at fixed spin 



FIG. 2: A schematic representation of the Parisi-Franz po- 
tential m, H^ for a model in the RFOT class. 




through the Franz-Parisi potentiaL For the plaquette 
overlap, the Franz-Parisi potential Vj^'^^{q) is defined 
as in p^ . with qs{gi,a') replaced by 
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FIG. 3: Local plaquette-plaquette correlation functions 
4>p (t) (solid lines), and the corresponding "intra-state" 
thermodynamical overlap (dashed lines) computed with the 
cavity method [25:]. Different curves represent the same tem- 
peratures as in Fig. [T] 



variable made of three interacting spins, i.e. we consider: 

1 ^ 
^^"""W = ]^E(^'iW'^»2W'^*^.W'^^i(0)a,2(0)a,3(0)) , 
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(19) 

Because the average of a plaquette variable is 
{aiiai2ai3) = tanh(/3), it is better to consider the con- 
nected and normalized function 
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This correlation is shown in Fig. |3] for different tempera- 
tures; it displays a behavior slightly different from the one 
shown in Fig. [H but still below the dynamical transition 
T < Trf it does never relax and it reaches a plateau whose 
value can be computed through the cavity method [25| . 
similarly to what done in the previous case. This plateau 
value represents the overlap of plaquettes of spins be- 
tween configurations belonging to the same state. This 
behavior is not surprising in light of the understanding 
of the dynamical transition as the decomposition of the 
measure in many states, within the RFOT theory pic- 
ture. With such local dynamics, the system is not able 
to escape from the initial state and the cavity method 
allows us to access precisely this intra-state properties. 

Because the dynamical behavior of the spin and pla- 
quette overlap is qualitatively equivalent, one would like 
to understand whether, from a thermodynamic point of 
view, the plaquette overlap can be used as an order pa- 
rameter to constrain the system in one metastable state. 
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An explicit calculation of this quantity is possible, by ex- 
ploiting the mapping towards the non-interacting model 
discussed in Section Hi CI thanks to the crucial fact that 
the plaquette overlap in XORSAT becomes the spin over- 
lap in the KCM. Therefore ^/^'"(q) = V,^^^\q). The 
computation is done in Appendix |^ and it shows explic- 
itly that the function V^'^^{q) is a convex function at 

all temperatures, with a single minimum at q = tanh (/3) 
corresponding to the thermodynamic value of qp in the 
paramagnetic state. We conclude that V^'~'^{q) does 
not signal the presence of a metastable state and not 
even provide an interpretation for the change in the dy- 
namical behavior of the plaquette correlation in Eq. (J19p 
when T <Td- 

This example shows that in order to find a satisfactory 
thermodynamic description that accounts for the dynam- 
ical behavior it is important to select the correct order 
parameter which in the case of the XORSAT model is 
the spin overlap. Even if, a priori, imposing a given spin 
overlap between two configurations might not seem very 
different from fixing the plaquette overlap, this leads to a 
qualitative different behavior of the corresponding Franz- 
Parisi potential. The reason is readily found by thinking 
to the structure of the energy landscape of XORSAT, 
and in particular, as in the case under study, when there 
is a single ground state. Up to an overall constant the 
ground state has energy zero and above it there are N ex- 
cited configurations that violate exactly one interaction 
(i.e. one plaquette has value —1). These configurations 
are at extensive Hamming distance among each other 
and with the ground state, therefore their spin overlap 
with the ground state is very close to zero. However, 
in terms of plaquettes, they are at distance 1 from the 
ground state and at distance 2 among each other, there- 
fore their plaquette overlap with the ground state and 
among each other is very large. This implies that the 
requirement to fix a high value of the plaquette overlap 
in the Franz-Parisi potential does not select a state and 
on the contrary allows to sum over configurations be- 
longing to different states. For this reason the overlap 
qp{g_,s!) docs not work as an order parameter that sig- 
nals the emergence of many metastable states that trap 
the dynamics. This was also recognized in the context of 
finite-dimensional plaquette models |12| . 

The situation is somehow analogous to what happens 
if one considers a ferromagnetic Ising model and adds an 
infinitesimal field which couples to SiSj . Clearly, this field 
does not break the symmetry and is therefore unable to 
detect the phase transition by constraining the system in 
one of the two symmetry-breaking pure states. However, 
the dynamical correlations of SiSj would clearly show a 



plateau in the low temperature phase if the system is 
constrained into one of its pure states. 



B. Kinetically constrained model 

We turn now to the discussion of the kinetically con- 
strained model, introduced in Section IIIBl Due to the 
mapping onto XORSAT, non-trivial properties in the dy- 
namics of the KCM are expected to appear in spite of the 
simplicity of the energy function Eq. ([5]) . 

The equilibrium local spin-spin correlations for this 
model read: 

N 



Cf^'''W-^E(^a(i).a(0)) 






KCM) 



gf^^^ft)-tanh^(/3) 
CifCM(o)_tanh2(/3) 



(22) 



where {s\GSkcm) is the ground state of Hrcm defined 
in Eq. (ITi|) . Thanks to the duality transformation the 
ground state of Hrcm is mapped into that of HxOR and 
the time evolution operator as well. Moreover the observ- 
able s'^ is mapped in a^i a^2 ^^3 in the XORSAT problem. 
Thus, the spin correlation in the KCM, Eq. p2p is equiv- 
alent to the plaquette correlation in the XORSAT prob- 
lem, Eq. (dll): Cf^^'^it) = C^°^(t). This equivalence is 
confirmed numerically in Fig. |4] at two different temper- 
atures, T = 0.45 and T = 0.6. The correlation functions 
of the KCM have been obtained using Metropolis dynam- 
ics in discrete time for a system with N — 60000 spins. 
In the case of the KCM, the problem of generating the 
initial equilibrium configuration is of course completely 
trivial. 

Thus, the dynamics of the spins in the KCM is the 
same shown in Fig. |3] for the plaquettes in XORSAT. It 
reveals non-trivial features and in particular a dynam- 
ical transition at the same critical temperature T^ = 
0.5098(2). Such a transition cannot be understood in 
terms of a simple thermodynamical reasoning on the ba- 
sis of the Franz-Parisi potential for the spin overlap of 
the KCM. Indeed, as we already discussed, K^'"*^(g) — 
yXOR^q-^ and both are trivial functions of q, see Ap- 
pendix [X] Indeed, due to the non-interacting nature of 
the spins of the KCM, from the static point of view the 
overlap qs{s,s') = jj Sa=i ^as'a is not expected to select 
any state and it cannot work as a static order parameter. 
This phenomenology is very reminiscent of other kineti- 
cally constrained models. At first sight thus, in the KCM 
one is not able to understand the dynamic transition as 
the consequence of a complex thermodynamic behavior. 

However, the mapping in Eq. (J12p teaches us that there 
is a static order parameter that does this work, i.e. what 
corresponds to the overlap between spin configurations 




100000 



FIG. 4: Comparison between the (connected and normalized) 
spin-spin correlation function of the KCM, (f>^ ' (t) (points), 
with the plaquette-plaquette correlation (j>p'-^^{t) of XORSAT 
(solid lines). The grey and the pink curves refer to diflent 
temperatures, respectively T — 0.6 and T = 0.45. 



of the XORSAT model. Expressed in terms of the spins 



Sa, it becomes: 



N N 

qni{s,s:) = -Y,l[{s,s',)y 

i—l a—1 



(23) 



However this quantity is a highly non local observable 
and then it cannot be easily and naturally recognized, 
especially if the mapping is not explicitly given (note that 
in particular the yf depend on the particular instance of 
random XORSAT under consideration). 



IV. DISCUSSION 

A. Summary of the results 
for the Bethe lattice model 

In this work we have shown that it is possible to es- 
tablish a mapping between the equilibrium dynamics of a 
glassy mean-field model described by the Random First- 
Order Transition theory, i.e. the XORSAT model, and a 
Kinetically Constrained Model. 

We have discussed the role of the overlap function that 
is chosen in the framework of RFOT theory to work as 
order parameter and signal the presence of metastable 
states. Overlap functions that couple different observ- 
ables define different distances between configurations 
and thus might fail to select a particular state. This is 
the case of the plaquette overlap in the XORSAT model. 
This observable is directly mapped into the spin overlap 
in the KCM. Therefore, the thermodynamic properties, 
and in particular the Franz-Parisi potential, are trivial 
when the plaquette overlap is considered in XORSAT, or 
the spin overlap is considered in the KCM. 
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However the dynamics of all local correlations shows 
no relaxation below a critical temperature T^. Within 
the XORSAT model, this is signaled by a static compu- 
tation of the Franz-Parisi potential for the spin overlap, 
as usual in RFOT theory. For the KCM, thanks to the 
mapping back to the XORSAT model, one can recast the 
non-trivial dynamical behavior within a thermodynami- 
cal description whose basic degrees of freedom however 
are hidden non-local functions of the spins. 

Thus, this example illustrates that -at least in this 
case- one is able to resolve the apparent dichotomy that 
exist between static and dynamic properties tracing them 
back to the physics of RFOT theory. This is not obvious 
from the point of view of the KCM: if one were to study 
it without any knowledge on the mapping on the RFOT 
model, identifying the correct order parameter would be 
impossible. 



B. Finite dimensional models 

As we discussed in the introduction, the mapping we 
used for the Bethe lattice model was already used ex- 
tensively in a class of finite dimensional plaquette mod- 
els [a. [lll - [l3| . However, these plaquette models do not 
seem to be described by the RFOT theory, because the 
energy cost due to the presence of an interface between 
two low energy states does not grow with system size ]la - 

On the other hand, we have shown in section |ll] that 
the mapping between the p-spin model and the KCM 
holds quite generally, under the only requirements that 



i) each spin is involved in exactly p interactions, ii) each 
interaction involves exactly p spin, and in) the ground 
state is unique. Therefore, it might be possible to devise 
a finite-dimensional model (possibly in three or higher 
dimensions) that satisfies these requirements and at the 
same time is well described by the RFOT scenario. This 
could be possibly achieved by extensions of the so-called 
gonihedric models [32, 33] or of the three dimensional 
plaquette model considered in |34| . A numerical compu- 
tation of the Franz-Parisi potential for th e sp in overlap 
in these models, following the strategy of [33, ISg]' could 
be instructive. 

Such a model could be quite useful to investigate the 
relation between the RFOT and KCM pictures in a finite 
dimensional setting. We recall, however, that the map- 
ping implies that T^ = 0, because the free energy must 
be analytic at all temperatures, therefore a RFOT sce- 
nario with a positive Kauzmann temperature would not 
fit in this analysis. 

Whether a similar mapping can be constructed in more 
general models remains an open problem [22i] . Finally, it 
would be very interesting to check whether a KCM with a 
continuous transition [37] could be constructed, for which 
a similar mapping exists. 
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Appendix A: Franz-Parisi potential for independent spins 



Here we report the calculation of the Franz-Parisi potential for the plaquette overlap in XORSAT, Vj^'^^{q), that 
is mapped onto the same potential for the spin overlap in the KCM, Vg^'^^{q). Since in the KCM the spins are 
independent, the latter is trivially computed. We have: 



1 



V.^'^^q) = Vf^'Hq) = -^ E^"'"""""^^^ l°g [^ e-^^--(^')<5(, - ,.(.,.') 



(Al) 



vifhere the spin overlap is qs{^,s!) — jjTlii'^i'^'i ^^^ ^ system of free spins in magnetic field Ekcm{£) = ^'Yli'^i 
Exploiting the relation log x ~ lim„_j.o 9„a;" the Franz-Parisi potential becomes: 



K 



KCM 



(q) = - lim dn V 



.;3 2:.T, 



L r^ .... T^ a—1 i 



— — lim d„ — , 

«^o "/3(2cosh 



= — lim dn — , 

n^o /3(2cosh 



lim d. 



\ /dAi dA e-'^^« 2:^=1 Aa y^ 



/3(E, n+T.,,a O+f! !:,,„ A^^nrf 



)+l3i:^^aTT' 



N 



n^o "/3(2cosh/3)W 



^^^^^-PNqn\+N\oef{\) 



with 



/(A) = E.,.1 ....." e^(-+i:„r»)+0Ai:„rr» ^ ^^ ^Pr 3 cosh (/3(1 + At) 



^ef' 



2 cosh 



(/3(1 + A) 



.-;3 



2cosh /3(1- A) 



Developing for small n 

/(A)=2cosh(/3) 
Thus 



2 cosh ^ 



(e-'^log 2cosh(/3(l-A) 



s'^loa 



2 cosh ^(1 + A) ) 



1 



K (<?) = --? lim a„ / dAe 



/3 ri^O 



-l3Nqn\+T 



i 



e^l^ log 



2 cosh ;3(1-A) 



+ 6*- 



2 cosh ;3(1+A) 



N . 
— mm 
p X 



/3qX 



1 



2 cosh /3 



-^loe 



2 cosh ( 13(1 - X)] +e^log 2cosh f;3(l + A)] ) 



(A2) 



(A3) 



(A4) 



(A5) 



and 
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9(A) - j^^ (e^ tanh(/3(l + A)) - e~^ tanh(/?(l - A))) (A6) 

The function (|A5[) is a convex function at all temperatures, with a single minimum at q* = q{X = 0) = (tanh/3)^, 
such that V.^^^^^q*) = F{T) = -|;iog2cosh(^). 



